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Although universal continuous-variable quantum computation cannot be achieved via linear optics 
(including squeezing), homodyne detection and feed-forward, inclusion of ideal photon counting 
measurements overcomes this obstacle. These measurements are sometimes described by arrays of 
beamsplitters to distribute the photons across several modes. We show that such a scheme cannot be 
used to implement ideal photon counting and that such measurements necessarily involve nonlinear 
evolution. However, this requirement of nonlinearity can be moved 'offline', thereby permitting 
universal continuous-variable quantum computation with linear optics. 
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I. INTRODUCTION 

The chief attraction of quantum computation is the 
possibility of solving certain problems exponentially 
faster than any known method on a classical com- 
puter and a significant effort is underway to real- 
ize a physical quantum computer j^, |3|. Optical real- 
izations of a quantum computer are particularly appeal- 
ing because of the robust nature of quantum states of 
light against the effects of decoherence as well as the 
advanced techniques for state preparation, photon ma- 
nipulation, and photodetection. Both discrete-variable 
(qubit-based) |[ ^| and continuous-variable (CV) 
schemes offer significant potential as optical quantum 
computers. However, the lack of a strong optical nonlin- 
earity has been a considerable hurdle for optical quantum 
computation 

A proposal by Knill, Laflamme, and Milburn 
(KLM) || describes how measurement of photons can 
be employed to induce a nonlinear transformation in a 
qubit-based optical quantum computer and how this pro- 
cedure can be done efficiently in a non-deterministic way. 
These remarkable results suggest that measurement in a 
CV system may be used to induce nonlinear evolution as 
well (although measurements of CV observables may not 
be possible in the von Neumann sense; see ||). A scheme 
proposed by Gottesman, Kitaev, and Preskill (GKP) ||] 
uses photon number measurement to induce a nonlinear 
transformation. If photon counting overcomes the obsta- 
cle of creating optical nonlinearities, then CV quantum 
computation may be feasible, just as the qubit-based lin- 
ear optical quantum computer may be feasible [^). Here 
we present three key results: 

1. universal quantum computation over continuous 
variables can be achieved using linear optics, homo- 
dyne measurement with feed-forward, and photon 
counting, 



2. the desired photon counting projective measure- 
ment cannot be performed using linear optics and 
existing photodetectors, and necessarily involves an 
optical nonlinearity, and 

3. the nonlinear transformations can be brought 'of- 
fline' to prepare quantum resources for a linear op- 
tical CV quantum computation to succeed in a de- 
terministic way. 

This paper is outlined as follows. We first review the 
Clifford group for continuous variables, which consists of 
linear optics transformations (including squeezing). The 
construction of a cubic phase state of GKP is outlined as 
a possible means to implement a nonlinear transforma- 
tion, necessary for universal CV quantum computation. 
This gate requires photon counting measurements, and 
we demonstrate that linear optics, homodyne measure- 
ment, and ideal photodetection is insufficient to imple- 
ment such schemes. A nonlinear interaction is shown to 
be a necessary component of any photon counting mea- 
surement. An analysis of truncated Hilbert spaces is in- 
cluded, and the paper concludes with a discussion on the 
use of photon counting measurements (and their associ- 
ated nonlinear transformations) offline in a CV quantum 
computation. 



II. CLIFFORD GROUP TRANSFORMATIONS 

The requirement of nonlinear transformations for uni- 
versal CV quantum computation can be understood 
by considering n harmonic oscillators, corresponding to 
n independent optical field modes, with annihilation op- 
erators {df, i = 1, . . . , n}. Linear optical transformations 
of these modes are described by unitary phase-space 
displacements (by mixing with 'classical fields' at beam 
splitters) 



Di(a) = exp(aa| — a*di) . 



(1) 



2 



with a e C; these transformations comprise the 
Heisenberg-Weyl group HW(n). For a classical pump 
field, parametric amplification invokes one-mode squeez- 
ing operations 

5 l (r/)=exp(i(^ l t2 -^a 2 2 )), (2) 

and two-mode squeezing operations 

Sij(r)) = exp(|(?7atat - rfo-idj)) , (3) 

for rj 6 C ||. (Although squeezing utilizes an optical 
nonlinearity of order two or higher, the transformation 
is regarded as being linear because the resultant Heisen- 
berg operator equations of motion are linear.) Squeezing 
operations (both one- and two-mode) generate the sym- 
plectic group Sp(2n, R). 

The squeezing operation £,-(77) with 77 real maps the 
canonical position as 

Si(y) ■ qi = ^fhj2{a % + at) -> exp(-r?)% ; (4) 

thus, the infinitely squeezed displaced vacuum 

lim S t { V )D t (q/V2h)\0) , (5) 

rj — >oo 

with q £ R is the (unnormalizable) position eigenstate 
\q)i. These position eigenstates are often employed as a 
computational basis for CV quantum computation, and 
are approximated in experiment by finite squeezing 0. 
Two-mode squeezing Sij(rj) acts in a similar fashion on 
the normal and antinormal modes ± dj, and the in- 
finitely squeezed two-mode vacuum 

|e> y = lim S^M, (6) 

rj— >oo 

with 77 € M is the EPR state satisfying 

yteg'ie^- = %-</)• (7) 

Two-mode squeezing also allows us to implement a uni- 
tary SUM gate g, defined as 

SUMy = exp(-^qiPj) = exp(A(aj + a,)(a] - %)) . 

(8) 

This gate acts on the computational basis of position 
eigenstates according to 

SXJM tj : \qi)i\qj)j -> \qi)i\qi + q 3 )j ■ (9) 

The « th mode is referred to as the control and the j th 
mode as the target. 

Phase-space displacements and squeezing together 
close to a finite-dimensional group known as the Clifford 
group (|, For n modes, the Clifford group is the 

semidirect product group [Sp(2n, R)]HW(n) generated 
by all Hamiltonians that are inhomogcneous quadratics 
in the canonical operators {ai,a\,i = 1,... , n}. The 



above unitary representation of the Clifford group is a 
subgroup of all unitary transformations on n modes. As 
such, they are insufficient to generate arbitrary unitary 
transformations and thus cannot perform universal quan- 
tum computation. The addition of a nonlinear operation 
such as that provided by the x^ 3 \ or optical Kerr, nonlin- 
earity JiT| suffices, in principle, to perform universal CV 
quantum computation, but is not feasible in quantum 
optical implementations due to the lack of sufficiently 
strong nonlinear materials with low absorption. How- 
ever, as stressed by Lloyd and Braunstein Q, any non- 
linear coupling on a single mode could allow for universal 
CV quantum computation, as opposed to the qubit case 
where a nonlinear coupling between qubits is required. 



III. THE CUBIC PHASE GATE 

One is lead to ask whether measurements can be used 
to induce nonlinear evolution, following the related ex- 
ample for qubit-based linear optics quantum computa- 
tion ||. First, we consider Clifford group transforma- 
tions conditioned on the results of projective-valued mea- 
surements (PVMs) in the computational basis (i.e., von 
Neumann measurements in the basis {\q),q € R})- As 
shown in [Tc| ], such measurements and feed-forward are 
efficiently simulatable on a classical computer and thus 
(under the assumption that universal quantum compu- 
tation is not efficiently simulatable classically) are also 
insufficient for universal quantum computation. These 
results also include a more realistic computational basis 
of finitely-squeezed states and realistic homodyne mea- 
surement of quadratures. Thus, in the following, we will 
consider Clifford group transformations conditioned on 
homodyne measurement to be part of "linear optics" . 

Possibly, measurements in a different basis can be 
employed to induce a nonlinear transformation. Such 
a scheme has been proposed (in the context of quan- 
tum computation with finite-dimensional qudits rather 
than CV) by GKP using measurement of photon num- 
ber. Specifically, these measurements are described by a 
PVM 

{En = |n)<n|,n = 0,1,2,...}, (10) 

for a single oscillator, where \n) is the eigenstate of the 
number operator N = a' a with eigenvalue n. In what fol- 
lows, we refer to this PVM as the photon counting PVM. 
The scheme of GKP is briefly outlined in the following, 
and relies on the creation of a so-called "cubic phase 
state" I7), which is the (unnormalizable) state defined as 

It) = J d i ex p(*7? 3 )k) • (11) 

The cubic phase state I7) can be prepared using 
squeezing, phase space displacement, and photon count- 
ing. Consider the two-mode squeezed vacuum state 
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Sl2{T})\tyi V G and a large momentum displacement 
of the first mode, to obtain the two-mode state 

\w,r)) = D x (iw)Su(ri)\0) , (12) 

with w £ R. By performing a measurement of photon 
number (described by the photon counting PVM) on the 
first mode, a measurement result of n photons projects 
the second mode of the pair into a cubic phase state I7') 
to a good approximation if w is sufficiently large (details 
can be found in |j, ^|), where 7' oc n^ 1 / 2 . This state 
can be transformed into \y) with 7 of order unity using 
one-mode squeezing. 

The cubic phase state can be used to implement a non- 
linear transformation on an arbitrary state of an op- 
tical mode as follows §. A SUM" 1 gate (see Eq. (|)) 
is executed with \ip)i as the control and as the tar- 
get. A position measurement is performed on the target, 
projecting the control into the state 

\i>') i = j {q = a\$VM7^) i \ 1 ) j 

= cxp(i(q i + a) 3 )\if;) l , (13) 

for a measurement outcome a. Invoking the Clifford 
group transformation 

U(a) = exp(iqf - i(q. t + a) 3 ) 

= exp(-ia 3 /4) exp(-3ai(^ + a/2) 2 ) , (14) 

(which can be implemented using linear optics) on the 
state gives a net transformation equivalent to ap- 
plying 

V 1 = exp(z 7 gf ) , (15) 

on \ip)i- We refer to the transformation V 7 , implemented 
in this manner, as the cubic phase gate. This nonlin- 
ear transformation could be used in combination with 
Clifford group operations to perform universal quantum 
computation over continuous variables. 

What is fundamentally important about this result is 
that all transformations involved, including 

1. the preparation of the two-mode squeezed state, 

2. the transformations needed to prepare I7), 

3. the SUM" 1 gate, and 

4. the transformation U(a) (which is quadratic in g), 

are implementable using Clifford group (linear optics) 
transformations. The resulting cubic phase gate is fully 
deterministic. The key component that allows for the 
nonlinear transformation is the measurement of photon 
number. In the following, we argue that such a measure- 
ment possesses "hidden" nonlinear evolution (i.e., the 
equivalent of an optical Kerr or higher order nonlincar- 
ity), and we can make this nonlinearity explicit. 



IV. PHOTON COUNTING 

Due to the classical simulatability results of CO), any 
CV quantum information process that initiates with 
finitely- or infinitely-squeezed vacua and employs only 
Clifford group transformations (phase space displace- 
ments, and one- and two-mode squeezing), homodyne 
detection, and classical feedforward, can be simulated ef- 
ficiently on a classical computer. If the cubic phase gate 
leads to universal quantum computation with continu- 
ous variables, then it must not satisfy the conditions of 
this theorem. (Otherwise, quantum computation could 
be simulated efficiently on a classical computer, which is 
believed to be impossible.) The key component of the cu- 
bic phase gate that does not satisfy the conditions of this 
theorem is the photon counting PVM; thus, we conclude 
that the photon counting PVM cannot be implemented 
using only Clifford group transformations and homodyne 
measurement. 

In the following, we show that this transformation can- 
not be implemented even with the addition of photode- 
tectors. We demonstrate that the photon counting PVM 
necessarily requires nonlinear evolution. 

A. Photon counting using threshold detectors 

The photon counting PVM employed by GKP consists 
of projections in the Fock state basis {|n), n = 0, 1, . . .}. 
For such a measurement to be performed, one requires 
photodetectors that can measure the number of photons 
in a mode; i.e., distinguish the state \n) from \n'), n ^ n' . 
We refer to such a photodetector as discriminating. How- 
ever, such discriminating photodetectors do not yet ex- 
ist jl2| . All photodetectors in current use effectively mea- 
sure whether there are no photons (n = 0), or at least 
one photon (n > 0) in a mode. (Note that a discrimi- 
nating photodetector, on the other hand, must be able 
to count intrinsically indistinguishable photons, that is, 
propagating photons in the same temporal, spatial and 
polarization mode [^3| .) We refer to existing photodetec- 
tors as threshold detectors, and a unit-efficiency detector 
as an ideal threshold detector (ITD). The PVM for an 
ITD is 

{n = |o>(o|, n >0 = /-|o)(o|}, (ie) 

where / is the identity operator. The second projector 
n>o projects onto the infinitely large subspace of states 
with one or more photon, due to saturation of the thresh- 
old detector @. 

It is possible to use ITDs to distinguish up to a finite 
number k of photons by using linear optics to couple the 
input mode to multiple ITD modes Q. For example, 
one could use an array of beamsplitters to distribute the 
photons of the input state over N modes such that it 
is highly unlikely that more than one photon is in any 
of the N modes f| ITDs are then used at each 
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mode, and the probability of undercounting photons is 
at most k(k - 1)/2N. For small k (as in KLM), the 
photon counting PVM can be approximated with high 
probability by using a sufficiently large N. 

In a related fashion, the visible light photon counter 
(VLPC) Jl6| has been constructed to discriminate be- 
tween one and two photons with a high degree of confi- 
dence, but the measurement does not correspond to the 
photon counting PVM; rather the VLPC is effective at 
distributing photons throughout the photosensitive re- 
gion, with localized regions acting as threshold devices 
(as with standard photodetectors). In other words, the 
single-mode input field is distributed amongst many lo- 
calised modes in the photodetector, each region operating 
as a threshold photodetector. Consequently the VLPC 
has much in common with the proposed detection of pho- 
tons via arrays of beam splitters to split the signal field 
(as discussed above), with an ITD existing at each output 
port. 

Whereas the use of multiple ITDs and linear optics ap- 
proximates a discriminating photodetector if the Hilbert 
space can be truncated as in qubit-based linear optical 
quantum computation, this scheme breaks down for CV 
quantum computation. Without a priori knowledge of 
the maximum number of photons in a mode (for CV, this 
number is infinite), one would require an infinite number 
of auxiliary modes and ITDs. (Below, we discuss resource 
issues even if the Hilbert space is truncated.) Thus, the 
photon counting PVM cannot be performed using linear 
optics (Clifford group transformations), homodyne mea- 
surement, and a finite number of ITDs. 

B. Photon counting with nonlinear evolution 

In order to implement the photon counting PVM, it is 
illustrative to employ a model of photodetection based 
on homodyne measurement and nonlinear optics (i.e., 
by employing a Hamiltonian that is cubic or higher in 
the photon creation and annihilation operators). In this 
model, a measurement of the phase shift in a probe field 
allows for a quantum non-demolition measurement of 
photon number in the signal field |l7|] . Consider the 
probe field to be a coherent state with large amplitude. 
We interact this probe field with an arbitrary signal field 
in a Kerr medium with interaction Hamiltonian 

H\nt = hxN sigaal Np Iohe , (17) 

where \ is proportional to the third-order nonlinear 
susceptibility. After an interaction time i, homodyne 
measurement is then used to infer the phase shift 4> in 
the probe field. Although such a homodyne measure- 
ment does not project the probe field into a phase state, 
one can in principle project to a state with arbitrarily 
small uncertainty A(j) in phase. A particular value of 
(j) can be used to infer the photon number n — 4>/{xt) 
to the nearest integer value, with corresponding uncer- 
tainty An = Acf>/(xt)- However, the photon number is 



only obtained modulo N — 27r/(xt); the periodicity of 
the phase does not give a true photon number measure- 
ment H[. This measurement projects the signal field 
| ip) into the subspace spanned by the number states \rij), 
where nj = {4> + 2irj)/x- 

To implement the photon counting PVM without is- 
sues of periodicity, we can couple the signal field to a 
pointer with an unbounded domain, such as the posi- 
tion q or momentum p of a probe. For example, the 
radiation pressure on a mirror is proportional to the flux 
of (monochromatic) photons that strike it, and a suitable 
coupling Hamiltonian would be |l8| 

Hint — AV s i gna ig pro be , (18) 

which is also nonlinear. For a probe field initially in 
the momentum eigenstate \p = 0), after an interaction 
time t a measurement of momentum p of the probe col- 
lapses the probe field into a momentum eigenstate \p) 
and thereby the signal field into a number state with 
n = p/(Xt) (again to nearest integer value). 

The resulting photon number measurement in either 
scheme will carry with it an error (related to measure- 
ment precision, and converting from continuous to dis- 
crete quantities). GKP require that An <§; n 1 / 3 for a 
functioning cubic phase gate; this condition places limits 
on the acceptable measurement errors. 

Thus, measurement of photon number can be de- 
scribed as a nonlinear interaction plus homodyne mea- 
surement. This result gives insight into the reason why 
such measurements can induce a nonlinear transforma- 
tion. Specifically, this model of photon number measure- 
ment is excluded from the conditions for efficient classical 
simulation fT(| , 

C. Truncation of photon number 

Naturally, one should be suspicious of the periodicity 
involved in the first measurement scheme (employing the 
Kerr interaction of Eq. (|l7])) and also of the unbounded 
nature of the second scheme (employing Eq. ([l^)). Any 
physical realization of a CV quantum information process 
must have finite energy, and thus the Hilbert space can 
effectively be truncated at some highest energy with pho- 
ton number n max ■ The issue of energy arises in the second 
scheme as well, where the momentum increases without 
bound and the energy is not bounded above or below. 
Even before infinite energy becomes an issue, a relativis- 
tic description must be applied. Coupling to momentum 
instead of the position does not eliminate these difficul- 
ties: the displacement itself must also be bounded by the 
physical boundaries of the laboratory. Truncation of the 
Hilbert space is an option that, however, presents chal- 
lenges Q ; for example, the system is no longer described 
by continuous variables but rather by large-dimensional 
qudits. 

An advantage of a truncated Hilbert space is that the 
linear optics scheme involving a finite N ITD modes be- 
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comes well-defined. Implementing these photon count- 
ing measurement schemes thus becomes a matter of re- 
sources. Consider the linear optics measurement scheme 
for maximum photon number n max . The probability of 
undercounting k photons is at most k(k — 1)/2N for 
TV the number of ITD modes. Thus, for a fixed prob- 
ability of undercounting, the required number of ITD 
modes scales as N oc nf nax , although detector inefficien- 
cies somewhat complicate the issue. For the (nonlinear) 
photon number measurement scheme involving the Kerr 
interaction of Eq. (|j"7|), however, one does not require ad- 
ditional modes or detectors so this quadratic scaling does 
not apply. All that is needed for the Kerr interaction 
scheme is an increase in phase resolution that behaves as 
Acj) oc n"^. Similar resolution arguments apply to the 
position-pointer scheme. Thus, even if one were to con- 
sider a truncated Hilbert space, it may be more practical 
to employ a nonlinear measurement scheme rather than 
a multimode ITD scheme; this result may be true for the 
KLM and GKP schemes as well. For true CV quantum 
computation, however, multiple ITD arrays cannot suf- 
fice, and nonlinear evolution is a necessary component of 
photon number measurement. 

V. CONCLUSIONS 

Despite the need for nonlinear evolution for photon 
counting, the cubic phase gate has a considerable ad- 
vantage over the use of nonlinear interactions directly. 
Specifically, this gate can be used to remove the use 
of nonlinear materials from the computation and utilize 
them only in the preparation of cubic phase states. In 
other words, the photon number measurement can be 
performed "offline", and the cubic phase states can be 
viewed as a quantum resource to be prepared prior to the 
computation. This way, the states used in the the com- 
putation need not pass through any optical Kerr nonlin- 
earities with their high absorption, thus avoiding the loss 
associated with using such materials. Also, if the proce- 
dure for producing cubic phase states possesses noise or 
other sources of error, imperfect cubic phase states can 



be purified to produce a smaller number of states with 
higher fidelity || . Again, an advantage is that this purifi- 
cation can be done offline and is not part of the computa- 
tion. This concept is similar to the KLM scheme, where 
the "difficult" gates are implemented offline on suitable 
ancilla states and then quantum teleported onto the en- 
coded states when needed. In our scheme, one simply 
prepares a sufficient number of cubic phase states prior 
to the computation, and the entire process may then oc- 
cur using only linear optics and homodyne measurement. 
A key advantage of this scheme is that the teleportation 
can be performed deterministically. 

In summary, we have shown that universal CV quan- 
tum computation can be obtained using linear optics 
(phase-space displacements and squeezing), homodyne 
measurement with classical feed-forward, and a realiza- 
tion of the photon counting PVM. We describe the PVM 
for current (ideal) photodetectors, and demonstrate that 
such detectors cannot be used to implement the pho- 
ton counting PVM with linear optics alone. This photon 
counting PVM carries with it implicit nonlinear evolu- 
tion, and we discuss how it can be implemented in a CV 
system using a Kerr interaction (or another nonlinear 
Hamiltonian) and homodyne measurement. The resource 
requirements of this measurement scheme compared with 
using linear optics and current photodetectors are out- 
lined. Finally, an advantage of this scheme is its use in 
the nonlinear gate of GKP, which removes the nonlin- 
ear operations from the computation and reduces them 
to "offline" preparation of ancilla states. These results 
place the implementation of strong nonlinear CV quan- 
tum gates, and thus universal CV quantum computation, 
in the realm of experimental accessibility. 



Acknowledgments 

This project has been supported by an Australian Re- 
search Council Large Grant. SDB acknowledges the sup- 
port of a Macquarie University Research Fellowship. We 
thank E. Knill and Y. Yamamoto for useful discussions. 



[1] P. W. Shor, "Proceedings, 35* Annual Symposium on 
Foundations of Computer Science," IEEE Press, Los 
Alamitos, CA (1994). 

[2] A. Steane, Rept. Prog. Phys. 61, 117 (1998). 

[3] M. A. Nielsen and I. L. Chuang, "Quantum Computa- 
tion and Quantum Information," Cambridge University 
Press, Cambridge (2000). 

[4] I. L. Chuang and Y. Yamamoto, Phys. Rev. A 52, 3489 
(1995). 

[5] E. Knill, R. Laflamme, and G. J. Milburn, Nature 409, 
46 (2001). 

[6] D. Gottesman, A. Kitaev and J. Preskill, Phys. Rev. A 
64, 012310 (2001). 



[7] S. Lloyd and S. L. Braunstein, Phys. Rev. Lett. 82, 1784 
(1999). 

[8] M. Ozawa, in "Quantum Communication, Computing, 
and Measurement," eds. O. Hirota, A. S. Holevo and C. 
M. Caves, Plenum Press, New York (1997), p. 233. 
[9] C. M. Caves, C. Zhu, G. J. Milburn, and W. Schleich, 
Phys. Rev. A 43, 3854 (1991). 
[10] S. D. Bartlett, B. C. Sanders, S. L. Braunstein and 

K. Nemoto, Phys. Rev. Lett. 88, 097904 (2002). 
[11] G. J. Milburn and D. F. Walls, Phys. Rev. A 28, 2065 
(1983). 

[12] K. J. Resch, J. S. Lundeen, and A. M. Steinberg, Phys. 
Rev. A 63, 020102(R) (2001). 



6 



[13] K. J. Blow, R. Loudon, S. J. D. Phoenix, and T. J. Shep- 
herd, Phys. Rev. A 42, 4102 (1990). 

[14] P. L. Kelley and W. H. Kleiner, Phys. Rev. 136, A316 
(1964). 

[15] P. Kok and S. L. Braunstein, Phys. Rev. A 63, 033812 
(2001). 

[16] J. Kim, S. Takeuchi, Y. Yamamoto, and H. H. Hogue, 
Appl. Phys. Lett. 74, 902 (1999). 



[17] B. C. Sanders and G. J. Milburn, Phys. Rev. A 45, 1919 
(1992). 

[18] A. F. Pace, M. J. Collett, and D. F. Walls, Phys. Rev. A 

47, 3173 (1993). 
[19] V. Buzek, A. D. Wilson-Gordon, P. L. Knight, and W. 

K. Lai, Phys. Rev. A 45, 8079 (1992). 



